A gauge invariant dressed holon and spinon description of the normal-state of 

underdoped cuprates 
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A partial charge-spin separation fermion-spin theory is developed to study the normal-state prop- 
erties of the underdoped cuprates. In this approach, the physical electron is decoupled as a gauge 
invariant dressed holon and spinon, with the dressed holon behaving like a spinful fermion, repre- 
senting the charge degree of freedom together with the phase part of the spin degree of freedom, 
while the dressed spinon is a hard-core boson, representing the amplitude part of the spin degree of 
freedom. The electron local constraint for the single occupancy is satisfied. Within this approach, 
the charge and spin dynamics of the underdoped cuprates are studied based on the t-t'-J model. It 
is shown that the charge dynamics is mainly governed by the scattering from the dressed holons due 
to the dressed spinon fluctuation, while the scattering from the dressed spinons due to the dressed 
holon fluctuation dominates the spin dynamics. 
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I. INTRODUCTION 

After over fifteen years of intense experimental stud- 
ies of doped cuprate superconductors, a significant body 
of reliable and reproducible data has been accumulated 
by using many probes 1 ' 2 , which shows that the most re- 
markable expression of the nonconventional physics is 
found in the normal-state 1,2 . The normal-state prop- 
erties exhibit a number of anomalous properties in the 
sense that they do not fit in the conventional Fermi- 
liquid (FL) theory 2-4 , and are closely related to the fact 
that these materials are doped Mott insulators 1-4 . The 
ground state in the undoped case is an antiferromagnetic 
(AF) long-range-order (AFLRO) Neel state, but chang- 
ing the carrier concentration by ionic substitution or in- 
creasing the oxygen content turns these compounds into 
strongly correlated metals leaving the AF short-range- 
order correlation still intact 1,2 . It is then not surpris- 
ing that the nonconventional behaviors are most strik- 
ing in the underdoped regime, where the concentration 
of doped holes is small. The anomalous properties ob- 
served in a variety of experiments, such as the nuclear 
magnetic resonance (NMR), nuclear quadrupole reso- 
nance (NQR), muon spin rotation (fiSR) techniques, in- 
elastic neutron scattering studies 5-10 , optical and trans- 
port measurements 11-14 , and angle-resolved photoemis- 
sion spectroscopy (ARPES) investigation 15 , exclude con- 
ventional theories. The single common feature in doped 
cuprates is the two-dimensional (2D) CuC>2 plane 1,2 , and 
it seems evident that the nonconventional behaviors are 
dominated by this plane. Very soon after the discovery 
of high-temperature superconductivity (HTSC) in doped 
cuprates, Anderson proposed a scenario of HTSC based 
on the charge-spin separation (CSS) in 2D 3 , where the 
internal degrees of freedom of the electron are decoupled 
as the charge and spin degrees of freedom, while the ele- 
mentary excitations are not quasi-particles but collective 
modes for the charge and spin degrees of freedom, i.e., 
the holon and spinon, then these holon and spinon might 
be responsible for the nonconventional behaviors. Many 



unusual properties of the underdoped cuprates are exten- 
sively studied following this line within the 2D t-J type 
model 4,16 . 

The decoupling of the charge and spin degrees of 
freedom of electron is undoubtedly correct in one- 
dimensional (ID) interacting electron systems 17 , where 
the charge and spin degrees of freedom are represented by 
boson operators that describe the excitations of charge- 
density wave and spin-density wave, respectively. In par- 
ticular, the typical behavior of the non-FL, i.e., the ab- 
sence of the quasiparticle propagation and CSS, has been 
demonstrated theoretically within the ID t-J model 18 . 
Moreover, the holon and spinon as the real elementary 
excitations in ID cuprates have been observed directly 
by the ARPES experiment 19 . Therefore both theoreti- 
cal and experimental studies indicate that the existence 
of the real holon and spinon is common in ID 20 . How- 
ever, the case in 2D is very complex since there are many 
competing degrees of freedom 1,2 . As a consequence, 
both experimental investigation and theoretical under- 
standing are extremely difficult 1-4 . Among the unusual 
properties of the underdoped cuprates, a hallmark is the 
charge transport 11-14 , where the conductivity shows a 
non-Drude behavior at low energies, and is carried by x 
holes, where x is the hole doping concentration, while the 
resistivity exhibits a linear temperature behavior over a 
wide range of temperatures. This is a strong experimen- 
tal evidence supporting the notion of CSS, since not even 
conventional electron-electron scattering would show the 
striking linear rise of scattering rate above the Debye fre- 
quency, and if there is no CSS, the phonons should affect 
these properties 21 . Moreover, it has been argued that 
the most plausible source of the absence of phonon scat- 
tering and of pair- breaking effects seems to be CSS 21 , 
and further, a compelling evidence for CSS in doped 
cuprates has been found from the experimental test of 
the Wiedemann-Franz law, where a clear departure from 
the universal Wiedemann-Franz law for the typical FL 
behavior is observed 22 . On the numerical study front, the 
crossover from the FL to non-FL behavior with decreas- 
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ing the hole doping concentration near the Mott tran- 
sition has been found within the 2D t-J model 23 . Fur- 
thermore, it has been shown within the 2D t-t'-J model 
by using the exact diagonalization method that there is 
a tendency of holes to generate nontrivial spin environ- 
ments, this effect leads to a decoupling of the spin from 
charge 24 . In this formal theory with the gauge 

invariant holon and spinon, i.e., the issue of whether the 
holon and spinon arc real, is centrally important 25 . In 
this paper, we propose a partial CSS fermion-spin theory, 
and show that if the local single occupancy constraint is 
treated properly, then the physical electron can be de- 
coupled completely by introducing the dressed holon and 
spinon. These dressed holon and spinon are gauge in- 
variant, i.e., they are real in 2D. As an application of 
this theory, we discuss the charge and spin dynamics of 
the undcrdoped cuprates within the t-t'-J model, and 
the results are qualitatively similar to that seen in exper- 
iments. 

The paper is organized as follows. The framework of 
the partial CSS fermion-spin theory is presented in Sec. 
II. Within this theory, the single-particle dressed holon 
and spinon Green's functions of the t-t'-J model are cal- 
culated in Sec. Ill by considering the dressed holon- 
spinon interaction, where the dressed holon and spinon 
self-energies are obtained by using the equation of motion 
method. In Sec. IV, we discuss the charge transport of 
the underdoped cuprates. The incommensurate (IC) spin 
response of the underdoped cuprates is studied in Sec. V. 
Sec. VI is devoted to a summary and discussions. 



II. GAUGE INVARIANT DRESSED HOLON AND 
SPINON 

We start from the t-t'-J model defined on a square 
lattice as 3,26 , 



(2) 



h = -tY^cic^ + t'Y^cla 

ifjcr if a 



i+f <7 



(i) 



where fj = ±x, ±y, t = ±x ± y, C\ a (Cj CT ) is the electron 
creation (annihilation) operator, Si = CjaCi/2 is spin 
operator with a = (a x , a y , a z ) as Pauli matrices, and /i 
is the chemical potential. The t-t'-J model (1) is subject 
to an important local constraint C| CT Cj CT < 1 that a 
given site can not be occupied by more than one electron. 
In the t-t'-J model, the strong electron correlation mani- 
fests itself by this constraint, and therefore the crucial re- 
quirement is to impose this constraint 3,26 " 28 . It has been 
shown that this constraint can be treated properly within 
the fermion-spin theory 29 . In this section, we show that 
the dressed holon and spinon in the fermion-spin theory 
are gauge invariant. Following the discussions in [29], we 
decouple the constrained electron operator as, 



with the constraint a la a ia = 1> where the spinless 
fermion operator hi keeps track of the charge degree of 
freedom, while the boson operator Oi CT keeps track of the 
spin degree of freedom, then the Hamiltonian (1) can be 
rewritten as, 

H = -t^2 h ia \ a h\ +fl a l+flry + t'Yl h i a \o h \+r a i+™ 

ifja if a 

- V^2 h l h i + J ^2(hih\)Si ■ S i+fl {h l+fl h\ +fi ), (3) 

i if/ 

where the pseudospin operator Sj = a\aai/2. In this 
case, the electron constraint ^2 a C\ a Ci a = 1 — h\hi < 1 

is exactly satisfied, with n^ 1 = h\hi is the holon num- 
ber at site i, equal to 1 or 0. This decoupling scheme 
is called as the CP 1 representation 30 . The advantage of 
this formalism is that the charge and spin degrees of free- 
dom of the electron may be separated at the mean field 
(MF) level, where the elementary charge and spin exci- 
tations are called the holon and spinon, respectively. We 
call such holon and spinon as the bare holon and spinon, 
respectively, since an extra U(l) gauge degree of freedom 
related with the constraint a\ a ai a — 1 appears, i.e., 
the CP 1 representation is invariant under a local U(l) 
gauge transformation, 



hi 



(4) 



and then all physical quantities should be invariant with 
respect to this transformation. Thus both bare holon 
hi and spinon ai a are not gauge invariant, and they are 
strongly coupled by these U(l) gauge field fluctuations. 
In other words, these bare holon and spinon are not real. 

However, the constrained CP 1 boson a ia can be 
mapped exactly onto the pseudospin representation de- 
fined with an additional phase factor, this is because 
the empty and doubly occupied spin states have been 
ruled out due to the constraint a\^a^ + a^a^ = 1, 
and only the spin-up and spin-down singly occupied spin 
states are allowed, therefore the original four-dimensional 
representation space is reduced to a 2D space. Due 
to the symmetry of the spin-up and spin-down states, 



occupied) | 



and | empty) | 



are singly- 



occupied and empty spin- up, while | occupied) j = f 



1 



and | empty) j. 



are singly-occupied and empty 



spin-down states, respectively, thus the constrained CP 1 
boson operators dj CT can be represented in this reduced 
2D space as, 



a j = e l T | occupied) | j (occupied 
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a j. = e J$i | occupied) | i (occupied | 

= e^ (° j) =e^S+, (5b) 

where S~ is the S z lowering operator, while S + is the 
S z raising operator, then the constraint J2 a a ia a i<? — 
S^S^ + 5j + = 1 is exactly satisfied. Obviously, the 
bare spinon contains both phase and amplitude parts, 
and the phase part is described by the phase factor e % ^ ia , 
while the amplitude part is described by the spin oper- 
ator Si. In this case, the electron decoupling form (2) 
with the constraint can be expressed as, 

G A = hje^Sr, C a = h\e^S+, (6) 

with the local U(l) gauge transformation (4) is rewritten 

as, 

hi^hie ie \ -> $ fa + fc. (7) 

Moreover, the phase factor of the bare spinon e l * iCT can 
be incorporated into the bare holon, thus we obtain a 
new fermion-spin transformation from Eq. (6) as, 

Cn = h\^S7, C U = 4S+ (8) 

with the spinful fermion operator h ia = e~ l ® i,r hi de- 
scribes the charge degree of freedom together with the 
phase part of the spin degree of freedom {dressed holon), 
while the spin operator Si describes the amplitude part 
of the spin degree of freedom {dressed spinon). This elec- 
tron decoupling form (8) is called as a partial CSS since 
only the amplitude part of the spin degree of freedom 
is separated from the electron operator. These dressed 
holon and spinon are invariant under the local U{\) gauge 
transformation (7), and therefore all physical quantities 
from the dressed holon and spinon also are invariant with 
respect to this gauge transformation. In this sense, the 
dressed holon and spinon are real. The dressed holon 
carries a spinon cloud (magnetic flux), and is a mag- 
netic dressing 24 . In other words, the dressed holon car- 
ries some spinon messages, i.e., it shares some effects of 
spinon configuration rearrangements due to the presence 
of the hole itself. We emphasize that the dressed holon 
h^ = e~ % ^ ia hi is the spinless fermion hi (bare holon) 
incorporated in the spinon cloud e~ 4 * iCT (magnetic flux). 
Although in the common sense hi„ is not an real spin- 
ful fermion operator, it behaves like a spinful fermion. 
In correspondence with these special physical properties, 
we find that h\ a h ia — h\e % ^ za e~ l ^ ia hi — h\h t , which 
guarantees that the electron constraint, ^2 a Cj a Ci a — 
SffH^Sr + Srhi^Sf = h^S+Sr + S rS+) = 

1 — h\hi < 1, is always satisfied in analytical cal- 
culations. Moreover the double spinful fermion occu- 
pancy, h\ a h\_ = e l9i 'h\h\e l9i -" = and h la h t - a = 
e~ l ® ia hihie~ % ^ i ~ a = 0, is ruled out automatically. Since 
the spinless fermion hi and spin operators and S~ 



obey the anticommutation relation and Pauli spin alge- 
bra, respectively, it is then easy to show that the spinful 
fermion hi„ also obeys the same anticommutation rela- 
tion as the spinless fermion hi. 

Although the choice of the CP 1 representation is con- 
venient, so long as h\hi = 1, J2 a Cj a Ci a = 0, no mat- 
ter what the values of S^S^ and 5r S^~ arc, therefore a 
"spin" even to an empty site has been assigned. It has 
been shown 29 that this defect can be cured by introduc- 
ing a projection operator Pi, i.e., the constrained elec- 
tron operator can be mapped exactly using the fermion- 
spin transformation (8) defined with an additional pro- 
jection operator Pi. However, this projection operator 
is cumbersome to handle in the many cases, and it has 
been dropped in the actual calculations 29 ' 31 . It has been 
shown 29,31,32 that such treatment leads to errors of the 
order x in counting the number of spin states, which is 
negligible for small dopings. Moreover, the electron con- 
straint is still exactly obeyed even in the MF approx- 
imation (MFA), and therefore the essential physics of 
the gauge invariant dressed holon and spinon is kept. 
This is because the constrained electron operator d a in 
the t-J type model can also be mapped onto the slave- 
fermion formalism 31 as Ci a — h\bi a with the constraint 
h\hi + J^cr "igbic = 1. We can solve this constraint by 
rewriting the boson operators bi a in terms of the CP 1 bo- 
son operators a ia as 6 icr = a icr \Jl — h\hi supplemented 

by the constraint X)o- a L a «'' = 1- As mentioned above, 
the CP 1 boson operators and with the constraint 
can be identified with the pseudospin lowering and rais- 
ing operators, respectively, defined with the additional 
phase factor, therefore the projection operator is ap- 
proximately related to the holon number operator by 

Pi ~ \J 1 — h\ a hi a = \J \ — h\hi, and its main role is 
to remove the spurious spin when there is a holon at the 
site i. 

III. DRESSED HOLON AND SPINON GREEN'S 
FUNCTIONS 

Before discussing the charge and spin dynamics, let 
us first calculate the dressed holon and spinon Green's 
functions. The low-energy behavior of the t-t'-J model 
in the partial CSS fermion-spin representation can be 
expressed as 29 ' 31 , 

h = -tJ2(hnS+hl +m S7 +n + hi^hl^) 

if) 

+ t> Y,(hiiS+hl +n S7 + . + h tl Srhl +n S+ +f ) 

if 

~ ^^2 h L h i" + Jcff^Sj • S i+fj . (9) 

icr if} 

where J c s = (1 — x) 2 J, and x = {h\ a hi a ) = {h\hi) is the 
hole doping concentration. As a consequence, the kinetic 
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part in the t-t'-J model has been expressed as the dressed 
holon-spinon interaction, which dominates the essential 
physics of the underdoped cuprates. The one-particle 
dressed holon and spinon two-time Green's functions are 
defined as, 

g a {% — j,t — t') = -i6(t - t')([h ia {t), h] a (t')\) 

= ((h itT (t);h] <7 (t'))), (10a) 

D(i .7, t t') = -i6{t - t')([S+(t), S7(t')]) 

= ((S+(t);Sr(t'))), (10b) 

respectively, where (. . .} is an average over the ensemble. 

A. Equation of motion 

Since the dressed spinon operators obey Pauli al- 
gebra, our goal is to evaluate the dressed holon and 
spinon Green's functions directly for the fermion and 
spin operators in terms of the equation of motion 
method. In the framework of the equation of mo- 
tion, the time- Fourier transform of the two-time Green's 
function G(lo) = ((A;A' f )) LU satisfies the equation 33 , 
w{(A;Al)) u = ([A, A*]) + (([A,H\;A^)) U . If we define 
the orthogonal operator L as, [A, H] = (A — iL with 
([L, A^}) = 0, the full Green's function can be expressed 
as, 

G(lu) = G {0) (lo) + ^G^(lu){{L;L^))^ \lo), (11) 

where <r = {[A, A^]), and the MF Green's function 
G (0)_1 (w) = (w - C)A- It has been shown 33 that if the 
self-energy £(w) is identified as the irreducible part of 
((L;£t)} w , the full Green's function (11) can be evalu- 
ated as, 



-<h^{SfST +f + SrS+ f ), 



(13c) 



G{w) = 



w-C-S(w)' 



(12) 



with E(w) = {{L-lS))Z r /<;. In the framework of the 
diagrammatic technique, corresponds to the contri- 
bution of irreducible diagrams. 



B. The mean-field theory 



Within MFA, the t-t'-J model (9) can be decoupled 



as, 



Hmfa =H t + Hj - 8Nt X i4>i + 8A^'x202, (13a) 
Ht = Xit h l+nv hi ° ~ ^ H h i+f <? hi ° 



IT}<7 



(13b) 



Hj - Jeff ^2[-e(Sf S i+ - + S t Sf + ~) + SfSf +fj ] 



where the dressed holon's particle-hole parameters <j)\ — 
{h\ a hi + fj a ) and 4> 2 — (h\ a h i+ f a ), the dressed spinon cor- 
relation functions Xl = (^t ^i+rt) an< ^ ^2 = ^i+f )' 
and e = 1 + 2t<pi/J e g. Since AFLRO in the undoped 
cuprates is destroyed 34 by hole doping of the order ~ 
0.024, there is therefore no AFLRO in the doped regime 
x > 0.025, i.e., (Sf) — 0, and a disordered spin liquid 
state emerges. It has been argued that this spin liq- 
uid state may play a crucial role in the mechanism for 
HTSC 3,4 . In this paper, we focus on the normal-state 
properties in the doped regime without AFLRO. In this 
similar MF theory 35 of the t-J model based on 
the fermion-spin theory has been discussed within the 
Kondo-Yamaji decoupling scheme 36 , which is a stage one- 
step further than the Tyablikov's decoupling scheme 37 . 
In this MF theory 35 , the phase factor e l ^ ia describing 
the phase part of the spin degree of freedom was not 
considered. Following their discussions 35 , we obtain the 
MF dressed holon and spinon Green's functions in the 
present case as, 



gi 0) (k,^) = 

D<°>(k,w) = 



Bk 



(14a) 
(14b) 



respectively, where B k = Ai[2xf(e7k - 1) + Xi(7k - 
e)] - X 2 (2 X z 2 j' k - X2 ), Ai = 2ZJ eJf , A 2 = 4Zfat, 
7k = (l/Z) £ - e^, 7 ^ = (1/Z) Ef e Jk ' f , Z is the num- 
ber of the nearest neighbor or second-nearest neighbor 
sites, and the MF dressed holon and spinon excitation 
spectra are given by, 



6c = ^Xi7k - Zt'\2lL - M. ( 15a ) 
= A^fk) 2 + A 2 {i k ) 2 + A zlkl ' k 
+ A Al k + A bl ' k + A e , (15b) 



respectively, with A\ = aeX 2 (exi + Xi/2), A 2 = 
aX 2 2 x 2 , A 3 = -a\ 1 \ 2 {exl + e%2 + Xi/ 2 )> A i = 
-eX\ [a(xl + e X i/2) + (aCf + (1 - a)/(4Z) - ae X i/(2Z)) + 
(ad + (1 - a)/(2Z) - a X f/2)/2] + aAiA 2 (C 3 + 
e X2 )/2, A 5 = -3aX 2 X2 /(2Z) + aAiA 2 (x! + e X i/2 + 
CI), A 6 = X{[a.Ci + (1 - a)/(4Z) - ae X i/(2Z) + 
e 2 (aC 1 + (1 - a)/(2Z) - a X f/2)/2} + X 2 (aC 2 + 
(1 - a)/(2Z) - a X z 2 /2)/2) - aeAiA 2 C 3 , and the 
spinon correlation functions Xi = (&i > X 2 — 

(SfSf +f ), c 1 = (i/z 2 )j: fhV ,(st n s; + .), cf = 
(V^ 2 )£^<^^,), c 2 = (i/z 2 )^,(s+ +f s; +T ,), 

Ca = (l/^)Ef(^ + +^r+f>. and 

Cf = (l/Z)J2f( S i+f) S i+f)- In order not to violate the 
sum rule of the correlation function (S^~S^) = 1/2 in the 
case without AFLRO, the important decoupling param- 
eter a has been introduced in the MF calculation 35 ' 36 , 
which can be regarded as the vertex correction. All 
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the above MF order parameters, decoupling parameter 
a, and chemical potential /i are determined by the self- 
consistent calculation 35 . 



C. The dressed holon and spinon self-energies 

With the help of Eq. (12), the full dressed holon and 
spinon Green's functions of the t-t'-J model (9) are ex- 
pressed as, 



5<r(k,w) = 
D(k,w) = 



^-£ fe -4 2) (k,^) : 



(16a) 
(16b) 



respectively, where the second-order dressed holon 
self-energy from the dressed spinon pair bubble 
4 2) (k,w) = {{Lf\t);L^\t'))) w with the orthog- 
onal operator if ) = -t £\ h l+ ^(Sr + .S+ - Xi) + 
t' J2f hi+T<j{S~ +f Sf — X2), and can be evaluated as 31 , 



+ 



, I F,""(t,py) 

\ U! + LUp+p' — LUpi — £p+fc 
LO + LUpi — LUp+pi — £p+fc 
LO + LUpi + LUp+pi — £p+fc 

Fj h) (k,p,p') \ 

LO — LUp+pi — LOpi — £,p+k I 



+ 



(17) 



where 7 1 2 2 (k,p,p / ) = [(t7p< +p+k - ^7p'+ p +k) 2 + 
(*7p'-k - £'7p<-k) 2 ]' F [ h \k,p,p') = n F (£ p+ k)[n B (ujp,) - 
n B {uj p+p >)} + n B {iO p+p >)[\ + ns(wp')], F^ h) (k,p,p') = 
n-F {£p+k ) [riB (<*y +p ) - n B (uy )]+n B {u p > )[l + n B (uy + P )] , 
Fs h \k,p,p') = n F (£ p+k )[l + n B (Lo p+ p,) + n B (u> p >)] + 
n B {ujp,)n B {ujp + p>), F { i h) {k,p,p l ) = n F (£ p+k) [l + 
n B (uJp+ P ' ) + n B (uy )]-[l + n B (uy )}[l + n B (lo p+p > )] , and 
n B (ujp) and n F {£,p) are the boson and fermion distri- 
bution functions, respectively. This dressed holon self- 
energy is ascribed purely to the dressed holon-spinon in- 
teraction, and characterizes the competition between the 
kinetic energy and magnetic energy The calculation of 
the dressed spinon self-energy is quite tedious 38 , since 
our starting point is the dressed spinon MF solution 35 
within the Kondo-Yamaji decoupling scheme in subsec- 
tion B. From Eq. (11) and the MF dressed spinon 
Green's function (14b), the full dressed spinon Green's 
function satisfies the relation 38 , Lo 2 D(k,uj) = Bk + 
(([[S+(t),ff(i)],H(t)];S7(f))) M with [[St,H],H] k = 



LOk^ k ~ *r^, s ' ) . In the disordered liquid state without 
AFLRO, the dressed holon-spinon interaction should 
dominate the essential physics 38 . In this case, the orthog- 
onal operator for the dressed spinon can be selected 
from rj; s) as 38 , 

77, a 

+ h l+n+ai h i+fil ~ 2 0a)S,% +a 

+ [{2x1 + £Xi)Ai - X2M *a(fcJ T /h+aT 

a 

+ hl+uhil - a 

f,a 

+ h\ +f+ ~ ai K +n -2^)S+ +f+ ., (18) 

where a — fj,T, with t fj = t, (f>^ = <pi, and tf = —t', 
4>f = 4>2- Following [38], we obtain the dressed spinon 



self-energy £< 2) (k,o>) = ({Lf'(() ; £>'"(('))).,„, 

^»(K.-) = ^(l) 2 S7 f 2 ( k ,p,p')^ 



(s). 



-«t, 



PP 1 



F[ s \k,p, P ') 



\.LO + £p+p' — £p' — LO k + p 

F { 2 s \k lP ,p') \ 

w + £,p+p' ~ ip' + ^fc+p / 



(19) 



with F[ s \k,p,p') 

"F^p+p')! 1 - n F {^)\ - n B (uj k +p)[n F {£,p>) - n F (^ p+p ,)}, 
and F 2 (s) (fc,p,p') = n F (W)[! - M&O] + [1 + 
n B (Lj k+p )][n F (^p>) - n F (€p+ P ')]- Within the dia- 
grammatic technique, this dressed spinon self-energy 
si 2 ^(k, lu) corresponds to the contribution from the 
dressed holon pair bubble, and is consistent with our pre- 
vious result 38 . 



IV. CHARGE TRANSPORT 

Recently, the emergence and evolution of the metallic 
transport in doped cuprates have been extensively stud- 
ied by virtue of systematic transport measurements 14 . 
It is shown that the resistivity shows a crossover from 
the low temperature insulating-like to moderate temper- 
ature metallic-like behavior in the heavily underdoped 
regime (0.025 < x < 0.055), and a temperature linear 
dependence with deviations at low temperatures in the 
underdoped regime (0.055 < x < 0.15). These striking 
behaviors have been found to be intriguingly related to 
the AF correlation 14 . In this case, a natural question 
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is what is the physical origin of this transport transfor- 
mation from the insulating liquid state in the heavily 
underdoped regime to the unusual metallic state in the 
underdoped regime? In this section, we try to discuss 
this issue. Since the local constraint has been treated 
properly in the partial CSS fermion-spin theory, the ex- 
tra U(l) gauge degree of freedom related with the local 
constraint is incorporated into the dressed holon as men- 
tioned in Sec. II. In this case, the external electronic field 
only is coupled to the dressed holons, and the conductiv- 
ity is given by, 



ct(w) = - 



lmU h ((j) 



LO 



(20) 



with ILh(io) is the dressed holon current-current cor- 
relation function, and is defined as 11^ (i — t') = 
((jh(t)jh(t'))), where the current density of the dressed 
holons is obtained by taking the time derivation of 
the polarization operator with the use of the equa- 
tion of motion as 31 , j h = {exit/2) JY-^ f}hl +fj<T h ia - - 

(e*2*' 7^)E l fa^ h l+fa^- With the help of the full 
dressed holon Green's function (16a), the current-current 
correlation function is evaluated as, 



x g a {\t,iuj' m 



where 

iuj. n is the Matsubara frequency, 7^(k) 

2 



sin 



(21) 



k x (Xit - 



2x 2 i' cos k y y + shY k y (xit - 2x2*' cos fc x ) ]/ 4 - The full 
dressed holon Green's function can be expressed as fre- 
quency integrals in terms of the spectral representation, 



g a (k, iu> n ) 



duj A { «\\l,uj) 
2ir iuj n — lo 



(22) 



with the dressed holon spectral function (k, lo) — 
— 2Img (7 (k, lo). Substituting the Eqs. (22) and (21) into 
Eq. (20), and evaluating the frequency summation, we 
obtain the conductivity as 31 , 

^)=(?) 2 ^E^ 2 (k) r ^^(k^'+c) 

2 N J-oo 2tt 



ka 

xA^(^') nF{UJ ' +W) ~ nF{UJ,) - 



(23) 



We have performed a numerical calculation for (t(lo), and 
the results at x = 0.03 (solid line), x = 0.05 (dashed line), 
and x = 0.07 (dotted line) for t/J = 2.5 and t'/t = 0.15 
with T = are plotted in Fig. 1, hereinafter the charge 
e is set as the unit. For a comparison, the experimental 
result 11 of YBa2Cu307_a; is also plotted in Fig. 1 (in- 
set). Our results show that there is a low-energy peak 
at lo < 0.3i separated by a gap or pseudogap at 0.3t 
from a midinfrared band. This midinfrarcd band is dop- 
ing dependent, the component increases with increasing 



dopings for 0.3* < u> < l.Oi and is nearly independent of 
dopings for lo > l.Oi, however, the position of the mid- 
infrared peak is shifted to lower energies with increased 
dopings. This reflects an increase in the mobile carrier 
density, and indicates that the spectral weight of the 
midinfrared sideband is taken from the Drude absorp- 
tion, then the spectral weight from both low energy peak 
and midinfrared band represent the actual free-carrier 
density. For a better understanding of the optical prop- 
erties, we have made a series of calculations for a(ui) 
at different temperatures, and the results at x = 0.06 
with T = (solid line), T = 0.1J (dashed line), and 
T = 0.3J (dotted line) for t/J = 2.5 and t'/t = 0.15 
are plotted in Fig. 2 in comparison with the experi- 
mental data 11 taken from YBa2Cu307- x (inset). It is 
shown that <j(lo) is temperature-dependent for lo < l.Oi, 
and almost temperature- independent for lo > l.Oi. The 
peak at low energies broadens and decreases in the height 
with increasing temperatures, while the component in the 
low energy region increases with increasing temperatures, 
then there is a tendency towards the Drude-like behav- 
ior. The midinfrared band is severely suppressed with 
increasing temperatures, and vanishes at high tempera- 
tures, in qualitative agreement with experiments 11 ' 12 . 
Now we turn to discuss the resistivity, which is eval- 




co/t 

FIG. 1. The conductivity at x = 0.03 (solid line), x = 0.05 
(dashed line), and x = 0.07 (dotted line) with t/J = 2.5 
and t'/t = 0.15 in T = 0. Inset: the experimental result of 
YBa 2 Cu 3 07- :c taken from Ref. [11]. 




co/t 

FIG. 2. The conductivity at x = 0.06 in T = (solid line), 
T = 0.1J (dashed line), and T = 0.3J (dotted line) with 
t/J — 2.5 and t'/t = 0.15. Inset: the experimental result of 
YBa2Cu307_o, taken from Ref. [11]. 
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uated as, p = l/adc, with the dc conductivity ad c is 
obtained from Eq. (23) as adc = lim^^o cr(u>). This 
resistivity has been evaluated numerically, and the re- 
sults are plotted in Fig. 3 as a function of tempera- 
ture at x = 0.03 (solid line), x = 0.04 (dashed line), 
x = 0.05 (dotted line), and x = 0.06 (dash-dotted line) 
for t/J = 2.5 and t'/t = 0.15 in comparison with the 
experimental data 14 taken from La2- x Sr :E Cu04 (inset). 
Our results show obviously that the resistivity is char- 
acterized by a crossover from the moderate temperature 
metallic-like to low temperature insulating-like behavior 
in the heavily underdoped regime, and a temperature 
linear dependence with deviations at low temperatures 
in the underdoped regime. But even in the heavily un- 
derdoped regime, the resistivity exhibits the metallic-like 
behavior over a wide range of temperatures, which also 
is in qualitative agreement with experiments 14 . 

The perovskite parent compound of doped cuprates is a 
Mott insulator, when holes are doped into this insulator, 
there is a gain in the kinetic energy per hole proportional 
to t due to hopping, but at the same time, the spin cor- 
relation is destroyed, costing an energy of approximately 
J per site. Thus doped holes in a Mott insulator can be 
considered as a competition between the kinetic energy 
(xt) and magnetic energy (J). The magnetic energy J 
favors the magnetic order for spins and results in frus- 
tration of the kinetic energy, while the kinetic energy xt 
favors derealization of holes and tends to destroy the 
magnetic order. In the present partial CSS fermion-spin 
theory, the scattering of dressed holons dominates the 
charge transport, since the scattering rate is obtained 
from the dressed holon self-energy S^(k, w), while this 
self-energy is evaluated by considering the dressed holon- 
spinon interaction, and characterizes a competition be- 
tween the kinetic energy and magnetic energy. In this 
case, the striking behavior in the resistivity is intrigu- 
ingly related to this competition. In the heavily under- 
doped regime, the dressed holon kinetic energy is much 
smaller than the dressed spinon magnetic energy at lower 




FIG. 3. The electron resistivity as a function of temper- 
ature at x = 0.03 (solid line), x = 0.04 (dashed line), 
x = 0.05 (dotted line), and x = 0.06 (dash-dotted line) with 
t/J = 2.5 and t'/t — 0.15. Inset: the experimental result of 
La2-xSr 2 ,Cu04 taken from Ref. [14]. 



temperatures due to the strong AF correlation, where 
the dressed holons are localized, and the scattering rate 
from the dressed holon self-energy is severely reduced, 
this leads to the insulating-like behavior in the resistivity. 
With increasing temperatures, the dressed holon kinetic 
energy is increased, while the dressed spinon magnetic 
energy is decreased. In the region where the dressed 
holon kinetic energy is larger than the dressed spinon 
magnetic energy at moderate temperatures, the dressed 
holons can move in the background of the dressed spinon 
fluctuation, then the dressed holon scattering would give 
rise to the metallic-like behavior in the resistivity. Since 
the charge transport is governed by the dressed holon 
scattering, then x dressed holons are responsible for the 
electron conductivity. 



V. INCOMMENSURATE SPIN DYNAMICS 

The interplay between AF correlation and HTSC in 
doped cuprates is now well-established 1 , but its full un- 
derstanding is still a challenging issue. Experimentally, 
NMR, NQR, and /iSR techniques, particularly inelas- 
tic neutron scattering, can provide rather detailed infor- 
mation on the spin dynamics of doped cuprates 5-10 . It 
has been shown 7,8 that when AFLRO is suppressed, the 
IC magnetic correlation develops at a quartet of wave 
vector [n(l ± 5),ir] and [it, (1 ± 5)n], where the incom- 
mensurability 5 increases almost linearly with the hole 
doping concentration x at lower dopings, and saturates 
at higher dopings. These exotic features are fully con- 
firmed by the data both on the normal and supercon- 
ducting states 7,8 . It has been argued that the emergence 
of the IC magnetic correlation is due to dopings 39 . Al- 
though a sharp resonance peak at the commensurate AF 
wave vector has been observed in some optimally doped 
samples, this commensurate scattering is the main new 
feature that appears in the superconducting phase 40,8 . 
In this section, we only discuss the IC magnetic correla- 
tion in the normal-state. Within the present partial CSS 
fermion-spin theory, the spin fluctuation couples only to 
the dressed spinons, then the dynamical spin structure 
factor (DSSF) is obtained in terms of the full dressed 
spinon Green's function (16b) as 38 , 

5(k, u) = -2[1 + n B (Lu)]ImD(k, u) 
= -2[l + n B (w)] 

B fe Im4 2) (k,a;) 



[lu 2 -uj 2 k - RcE< 2) (k, uj)} 2 + [Im£^ (k, w)] 

(2) (2) 

where ImSs (k, uS) and ReSg (k, u) are the correspond- 
ing imaginary part and real part of the dressed spinon 
self-energy function Sg 2 ^(k, u>) in Eq. (19). 

At the half-filling, the spin fluctuation scattering re- 
mains commensurate at the AF wave vector Q = 
[1/2,1/2] position (hereafter we use units of [27r,27r]), 



(2), 



(24) 
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which is not presented here for the sake of space. In- 
stead, we plot the DSSF spectrum 5(k, to) in the (k x , k y ) 
plane at x = 0.06 with T = 0.05J and u = 0.05J for 
t/J = 2.5 and t'/t = 0.15 in Fig. 4. This result shows 
that with dopings, there is a commensurate-IC transition 
in the spin fluctuation geometry, where all IC peaks lie 
on a circle of radius of 5. Although some IC satellite di- 
agonal peaks appear, the main weight of the IC peaks is 
in the parallel direction, and these parallel peaks are lo- 
cated at [(l±S)/2, 1/2] and [1/2, (l±6)/2\. The IC peaks 
are very sharp at low temperatures and energies, which 
means that these low energy excitations have a dynam- 
ical coherence length at low temperatures that is larger 
than the instantaneous correlation length. For consider- 
ing IC magnetic fluctuation at a relatively high energy, 
we have made a series of scans for S(/k, ui) with several 
energies, and the result at x = 0.06 in t/J = 2.5 and 
t'/t = 0.15 with T = 0.05J for u = 0.1J is shown in 
Fig. 5. Comparing it with Fig. 4 for the same set of 
parameters except for lu = 0.05J, we see that at low 
temperatures, although the positions of the IC peaks are 
energy independent, the IC peaks broaden and weaken 
in amplitude as the energy increase, and vanish at high 
energies. This reflects that the excitation width increases 
with increasing energies, and thus leads to that the life- 
time of the excitations decreasing quickly with increasing 
energies. The present DSSF spectrum has been used to 
extract the doping dependence of the incommensurability 
S(x), which is defined as the deviation of the peak posi- 
tion from the AF wave vector position, and the result 
is plotted in Fig. 6 in comparison with the experimen- 
tal result 7 taken from La2-xSr x Cu04 (inset). Our result 
shows that S(x) increases progressively with the doping 




FIG. 4. The dynamical spin structure factor spectrum in 
the (k x ,ky) plane at x = 0.06 in T = 0.05J and uj = 0.05J 
with t/J = 2.5 and t'/t = 0.15. 




FIG. 5. The dynamical spin structure factor spectrum in 
the (k x , k y ) plane at x = 0.06 with t/J — 2.5 and t'/t = 0.15 
in T = 0.05J and lj = 0.1 J. 



concentration at lower dopings, but saturates at higher 
dopings, in qualitative agreement with experiments 7 ' 8 . 

The universal integrated dynamical spin response is 
a characteristic feature, and is closely related to many 
other normal-state properties 1 . The integrated dynami- 
cal spin response is manifested by the integrated dynam- 
ical spin susceptibility (IDSS), and is expressed as, 

^ r )4^ /(k ' w) ' (25) 

k 

where the dynamical spin susceptibility is related 
to DSSF by the fluctuation dissipation theorem as, 
X"(k,w) = (1 - e-^)S(k,uj) = -2Im£)(k, w). This 
IDSS has been evaluated numerically, and the results at 
x = 0.12 for t/J = 2.5 and t'/t = 0.15 with T = 0.2J 
(solid line), T = 0.3J (dashed line), and T = 0.5J (dot- 
ted line) are plotted in Fig. 7 in comparison with the 
experimental data 9 taken from La2- a; Sr :E Cu04 (inset). 
It is shown that the shape of IDSS appears to be par- 
ticularly universal, and can be scaled approximately as 
I(uj,T) (x arctan[aio;/T + a3(w/T) 3 ], where I(w,T) is al- 
most constant for lo/T > 1, and then begin to decrease 
with decreasing lo/T for lo/T < 1, also in qualitative 
agreement with experiments 9,10 . 

Although the scattering of the dressed spinons domi- 
nates the spin dynamics, the effect of the dressed holons 
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FIG. 6. The doping dependence of the incommensurability 
5(x). Inset: the experimental result of La2- a; Sr a ;Cu04 taken 
from Ref. [7]. 

0.4 i 1 




FIG. 7. The integrated dynamical spin susceptibility at 
x = 0.12 with t/J = 2.5 and t'/t = 0.15 in T — 0.2J (solid 
line), T = 0.3 J (dashed line), and T = 0.5 J (dotted line). 
Inset: the experimental result of L^-mSr^CuCU taken from 
Ref. [9]. 
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on the dressed spinon part is critical in determining 
the characteristic feature of the IC magnetic correla- 
tion, which can be understood from the properties of 
the dressed spinon excitation spectrum E% — u 2 + 

(2) 

ReEr'(k,£ fc ). During the calculation of DSSF spec- 
trum in Eq. (24), we find when W(kg,uj) = [oj 2 - w\ - 

, lo)] 2 ~ at some critical wave vectors ika in 
low energies, the IC peaks appear, then the weight of the 
IC peaks is dominated by the inverse of the imaginary 

(2) 

part of the dressed spinon self-energy 1/ImEs (ks,u>). 
Thus the positions of the IC peaks are determined by 

( 2) 

both functions W(k, uj) and IniSs (k, u>), where the zero 
points of W(k,uj) (then the critical wave vectors k^) is 
doping dependence. Near the half-filling, the zero point 
of W(k,u>) locates at the AF wave vector [1/2,1/2], so 
the commensurate AF peak appears there. With dop- 
ings, the holes disturb the AF background. Within the 
partial CSS framework, as a result of the self-consistent 
motion of the dressed holons and spinons, the IC mag- 
netic correlation is developed beyond a certain critical 
doping, this reflects the fact that the low energy spin 
excitations drift away from the AF wave vector, or the 
zero point of W(ks,u>) is shifted from [1/2,1/2] to kg. 
As seen from Eq. (24), the physics is dominated by 
the dressed spinon self-energy renormalization due to 
the dressed holon pair bubble. In this sense, the mo- 
bile dressed holons are the key factor leading to the IC 
magnetic correlation, i.e., the mechanism of the IC type 
of structure in doped cuprates is most likely related to 
the dressed holon motion. This is why the position of the 
IC peaks can be determined in the present study within 
the t-t'-J model, while the dressed spinon energy depen- 
dence is ascribed purely to the self-energy effects which 
arise from the dressed holon-spinon interaction. Since 

(2) 

the values of ImSs (k, w) increase with increasing ener- 

(2) 

gies, then all values of \/\mTr s (k, u>) are very small at 
high energies, which leads to the IC peaks disappearing 
at high energies. 

VI. SUMMARY AND DISCUSSIONS 

In summary, we have developed a partial CSS fcrmion- 
spin theory to study the physical properties of the under- 
doped cuprates. In this approach, the physical electron 
is decoupled completely as the dressed holon and spinon, 
with the dressed holon keeps track of the charge degree of 
freedom together with the phase part of the spin degree 
of freedom, while the dressed spinon keeps track of the 
amplitude part of the spin degree of freedom. The elec- 
tron local constraint for single occupancy is satisfied in 
analytical calculations. The dressed holon is a magnetic 
dressing, and it behaves like a spinful fcrmion, while the 
dressed spinon is neither boson nor fermion, but a hard- 
core boson. Moreover, both dressed holon and spinon are 
gauge invariant, and in this sense, they are real and can 



be interpreted as physical excitations. In the common 
decoupling scheme, we obtain the full dressed holon and 
spinon Green's functions by using the equation of motion 
method. 

Within this theoretical framework, we have studied the 
charge and spin dynamics of the undcrdoped cuprates 
based on the t-t'-J model. The conductivity spectrum 
contains a non-Drude low energy peak and a broad midin- 
frared band, while the temperature dependent resistivity 
is characterized by a crossover from the moderate tem- 
perature metallic-like to the low temperature insulating- 
like behavior in the heavily underdoped regime, and a 
temperature linear dependence with deviations at low 
temperatures in the underdoped regime. The commensu- 
rate neutron scattering peak at the half-filling is split into 
IC peaks with dopings, where the incommensurability is 
doping dependent, and increases with the hole doping 
concentration at lower dopings, and saturates at higher 
dopings. These results are qualitatively similar to those 
seen in experiments. It is essential that these theoreti- 
cal results were obtained without any adjustable param- 
eters. These results also show that the charge dynamics 
is mainly governed by the scattering from the dressed 
holons due to the dressed spinon fluctuation, while the 
scattering from the dressed spinons due to the dressed 
holon fluctuation dominates the spin dynamics. In this 
case, the spin and charge dynamics in the normal-state 
are almost independent, and the perturbations that in- 
teract primarily with the charge do not much affect the 
spin 21 , therefore the notion of partial CSS naturally ac- 
counts for the qualitative features of the normal-state of 
the undcrdoped cuprates. 

Based on this partial CSS fermion-spin theory, we have 
discussed the mechanism of HTSC in doped cuprates . 
It is shown that dressed holons interact occurring directly 
through the kinetic energy by exchanging the dressed 
spinon excitations, leading to a net attractive force be- 
tween the dressed holons, then the electron Cooper pairs 
originating from the dressed holon pairing state are due 
to the charge-spin recombination, and their condensation 
reveals the superconducting ground-state. The electron 
superconducting transition temperature is determined by 
the dressed holon pair transition temperature, and is pro- 
portional to the hole doping concentration in the under- 
doped regime. To our present understanding, the main 
reasons why the present theory is successful in study- 
ing the physical properties of doped cuprates are that 
(1) the electron local constraint is exactly obeyed during 
analytical calculations in contrast with the slave-particle 
approach, where the local constraint is explicitly replaced 
by a global one 2 ' 29 . In this case, the representation space 
in the slave-particle approach is much larger than the rep- 
resentation space for the physical electron. (2) Since the 
extra U(l) gauge degree of freedom related with the local 
constraint has been incorporated into the dressed holon, 
this leads to the fact that the dressed holon and spinon 
are gauge invariant in the partial CSS fermion-spin the- 
ory. However, the bare holon and spinon in the slave- 
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particle theory are strongly coupled by the U(l) gauge 
held fluctuation, they are not gauge invariant. (3) The 
representation of the dressed spinon in terms of the spin 
raising and lowering operators is essential in the present 
approach 29 , because whenever a dressed holon hops it 
gives rise immediately to a change of the spin background 
as a result of careful treatment of the constraint given in 
Sec. II. This is why the dressed holon-spinon interac- 
tion (kinetic part) dominates the essential physics of the 
undcrdoped cuprates 31 ' 38 . 
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